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Abstract 

The understanding of the kinetic properties of hydrogen (isotopes) adatoms on graphene is important 
in many fields. The kinetic properties of hydrogen-isotope (H, D and T) monomers were simulated using 
a composite method consisting of density functional theory, density functional perturbation theory and 
harmonic transition state theory. The kinetic changes of the magnetic property and the aromatic n bond of 



the hydrogenated graphene during the desorption and diffusion of the hydrogen monomer was discussed. 

The vibrational zero-point energy corrections in the activation energies were found to be significant, ranging 

from 0.072 to 0.205 eV. The results obtained from quantum-mechanically modified harmonic transition state 

£> ■ theory were compared with the ones obtained from classical-limit harmonic transition state theory over a 

ov 

wide temperature range. The phonon spectra of hydrogenated graphene were used to closely explain the 
(reversed) isotope effects in the pref actor, activation energy and jump frequency of the hydrogen monomer. 
The kinetic properties of the hydrogen-isotope monomers were simulated under conditions of annealing 



for 10 minutes and of heating at a constant rate (1.0 K/s). The isotope effect was observed; that is, a 
hydrogen monomer of lower mass is desorbed and diffuses more easily (with lower activation energies). 



The results presented herein are very similar to other reported experimental observations. This study of the 
kinetic properties of the hydrogen monomer and many other involved implicit mechanisms provides a better 
understanding of the interaction between hydrogen and graphene. 

PACS numbers: 68.65.Pq, 68.43.Bc, 68.35 Ja 



I. INTRODUCTION 

The fact that hydrogen interacts with graphite surface in outer space 1-4 and in fusion devices 4 - 5 
requires a clear understanding of the kinetic properties of hydrogen adatoms on graphite surfaces 
for astronomy and the nuclear industry. Such properties are also very important for the realization 
of both graphite-based hydrogen storage 6 and graphene-based electronics. -~— 

The desorption and diffusion of hydrogen adatoms on graphite surfaces are two important ki- 
netic processes whose properties depend sensitively on the interaction between the adatom and the 
graphite surface. Consequently, investigation of the desorption and diffusion of hydrogen adatoms 
on graphite surfaces can further our understanding of the mechanisms involved in these interac- 
tions. Some researches have been carried out on the kinetic properties of hydrogen adatoms on 
graphite surfaces.- 13 ' 16117 The adsorption of hydrogen monomers on graphene is an essential step in 
hydrogenation, and the kinetic properties of this step determine the outcome of the hydrogenation 
process. Hornekaer et. al. found that a large fraction of the D monomers on a graphite surface will 
be desorbed during annealing at room temperature for 10 minutes, while the D monomers on the 
graphite surface are diffusionally immobile at room temperature. 3 Such diffusional immobility is 
consistent with Baouche's observations in time-programmed desorption experiments. 18 Some the- 
oretical studies on the kinetic properties of hydrogen monomers have also been reported . 3i4il9 ~ — 
The calculated chemisorption energies of hydrogen monomer on graphene vary around a value of 
0.75 eV with a range of about 0.3 eV; this magnitude of variation represents quite a large ener- 
getic uncertainty in kinetics. Furthermore, in kinetic simulations, the vibrational characteristics 
of the light-mass hydrogen monomers should be considered. Thus, a fully theoretical simulation 
that can be precisely and directly compared with experimental observations in both desorption and 
diffusion is still lacking. 

In this report, the kinetic properties of various hydrogen-isotope (H, D and T) monomers were 
simulated by means of a composite method consisting of density functional theory (DFT), 23 den- 
sity functional perturbation theory (DFPT) 24 and harmonic transition state theory (hTST).— ~— 
The kinetic change of the magnetic property of the hydrogenated graphene during the desorption 
and diffusion of the hydrogen monomer was discussed. The vibrational contribution, including 
the zero-point energy correction in the activation energy, was considered with hTST The kinetic 
properties of various hydrogen-isotope monomers were simulated under conditions of annealing 
for 10 minutes and of heating at a constant rate (1.0 K/s). 



II. METHODOLOGY 

The over-barrier jump frequency between two local minimum states (initial and final states or 
reactant and product states) can be expressed in the Arrhenius form a s 27 ' 29 

E 
v = v*exp(-—*p) (1) 

k B T 

where v* is the exponential prefactor and E ac is the activation energy that is required for the reaction 
to occur. The activation energy is defined herein to be the vibrational zero-point energy corrected 
potential barrier, which is expressed as 

3AT-1 , 3N 



E ac = AE p + -J]^--Y J hco I i 

= AE p -AF vib (0) (2) 



where AE P is the potential barrier in the reaction path, which can be obtained from DFT calcu- 
lations, and to' and co^ are the vibrational frequencies of the z'fh mode in the initial and saddle- 
point states, respectively. Correction AF v a,(0) comes from the vibrational zero-point energy dif- 
ference between the initial and saddle-point states. The total vibrational degrees of freedom are 
3N; an imaginary vibrational mode along the migration coordinate in the saddle-point state is 
excluded; thus, 3N-1 vibrational modes are considered for the saddle-point state. From quantum- 
mechanically modified hTST, prefactor is expressed a s 25 ' 29 ' 30 
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where h T {a>i) is the bosonic phonon occupation number of the ith vibrational mode. When the 
temperature approaches infinity, the classical limit of the prefactor is expressed as 
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v '« = Tn^r- s (4) 

n <4 
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This classical-limit form is also the Vineyard's form,— in which all of the vibrational modes are 
assumed to be completely thermo-activated. From hTST as described above, the mass -dependent 
vibrational frequencies are responsible for the difference in kinetic properties of Ff, D and T 
monomers as shown below in this report. 



The first-order rate equation for desorption can be expressed as^^ 

^ = -v des (T)n(t) (5) 

at 

where v des is the desorption frequency; n(i) is the residual number of hydrogen monomers on 
graphene at the t moment, and t = represents the starting moment for initializing the kinetic 
movement. 

In the annealing process, the temperature is kept constant (T = To); in that case, the variation 
of the residual number with respect to time is expressed as 

n(t) = n(0) exp [-v des (T )t] (6) 

The diffusion can be described by the mean square displacement (< \r(t) - r(0)| 2 >) of a monomer 
parallel to graphene. From Fick's second law, we have 

<\r(t)-r(0)\ 2 >=2dD ad (T)t (7) 

where d is the dimensionality of the diffusion of a hydrogen monomer on graphene (taken to be 2 
here) and D ad is the diffusion coefficient, which is temperature-dependent and is expressed as 29 

D ad = ^a 2 (8) 

2d 

where T is the total jump frequency of the monomer and a is the jump length. For the diffusion 
of a monomer on graphene, Y is taken here to be 3v d iff (v d iff is the diffusion frequency), which is 
combined with the fact that there are three equivalent final sites to which a monomer at an initial 
site can diffuse, and a is taken to be the optimized C-C bond length of 1.426 A. The diffusion 
radius of the monomer is defined as the square root of the mean square displacement 



r diff = ^2dD ad (T)t (9) 

which can be used directly to determine whether the hydrogen monomer on graphene is diffusion- 
ally mobile. 

In the heating process, if the temperature increases at a constant rate a (T = at), the variation 
of the residual number with respect to time is expressed as 



Jo 



n(0 = n(0)exp[- v des (at)dt] (10) 



where a is always taken to be 1.0 K/s . 2il6il? The relative desorption rate is therefore defined by the 
ratio of the desorption rate to the initial monomer number 

R de s(T,t)=v des (TMt)/n(0) (11) 

which is both time- and temperature-dependent. 

In this report, a monolayer of graphene served as a structural model of a graphite surface. 
This was feasible because the effect of the weak van de Waals interaction between neighbor- 
ing graphene layers on the chemisorption properties of hydrogen monomer can be neglected; in 
our DFT tests using local density approximation (LDA) and generalized gradient approximation 
(GGA), the bilayer-graphene structural model yields potential barriers with a difference less than 
8 meV with respect to the monolayer structural model. The lattice constant of graphene has been 
found very close to graphite (within 0.05%) at 200 - 500 K, 31 which is the most important tem- 
perature range for the kinetic movements of the hydrogen monomer as shown below. A hydrogen 
monomer on a 5 x 5 periodic supercell of graphene (50 C atoms) with a 10 A vacuum along the 
direction perpendicular to the surface can be considered as an isolated monomer; our tests show 
that the interaction between two monomers in two such neighboring supercells can be neglected. 
This is consistent both with the findings of Shytov 7 based on a theory of electron-mediated in- 
teraction and with those of Boukhvalov 32 based on DFT calculations. The DFT and DFPT cal- 
culations were carried out using the Quantum Espresso code package, 33 in which the ultrasoft 34 
spin-polarized PBE 35 pseudopotentials were used to describe the electronic exchange and corre- 
lation energy. The wave function and the charge density were expanded using energy cutoffs of 
35 and 350 Ry. For the calculation of electronic density of states (DOS), a 6 x 6 x 1 uniform k- 
point grid with the tetrahedron interpolation scheme 36 was used for the integration of the electronic 
states over the first Brillouin zone. And for the other calculations, a 4 x 4 x 1 uniform &-point grid 
with the Methfessel-Paxton smearing technique 37 was used, where the smearing width was taken 
to be 0.03 Ry. The partial DOS of hydrogen was obtained using the Lowdin population analysis. 
A force threshold of 10" 4 Ry/bohr was used for the structural optimization. The reaction paths 
were described by the minimum energy paths (MEPs) between two local minimum states that 
were calculated using the climbing-image nudged elastic band method 38 with 5 images for each 
reaction path. For calculating the vibrational frequencies, only the gamma point at the Brillouin 
zone center was selected. 



III. RESULTS AND DISCUSSION 

The calculated MEPs for the desorption and diffusion processes of a hydrogen monomer are 
shown in Fig. \T\ in which the structures of the initial (reactant), saddle-point (transition) and final 
(product) states are also shown. The final state in the desorption MEP has been fixed to be the 
physisorption state. In order to simulate the desorption MEP from the physisorption state to the 
fully desorbed state, a vacuum with a size larger than 10 A (e.g. 16.6 A in Ref. |3|) is needed to 
completely eliminate any weak interaction between the desorbed hydrogen atom and the graphene 
layer in the neighboring supercell. However, as seen from the methodology part, precise potential 
surfaces from the chemisorption state to the saddle-point states are enough for the simulation of 
reaction rates in this report. The energy for the fully desorbed state here is calculated by summing 
the energies of an isolated hydrogen atom and an isolated graphene layer. In the desorption process 
shown in Fig. ED (a), the hydrogen monomer moves away from the graphene layer in the direction 
perpendicular to the layer. At the first stage of the desorption process (the chemisorption state), 
the hydrogen atom is at the top site with a height of 1.58 A bonded to a C atom that protrudes 
from the graphene plane with a height of 0.45 A. The calculated protrusion height is somewhat 
method-dependent and values ranging from 0.26 to 0.49 A have been obtained by others .-^i^- 
In the desorption saddle-point state, the hydrogen atom is at the top site with a height of 1.98 A. 
The potential barrier was defined as the energy difference between the saddle-point state and the 
initial state; the calculated desorption potential barrier is 1.075 eV. The final (physisorption) state 
is 0.865 eV higher in energy than the initial (chemisorption) state; the fully desorbed state is 0.890 
eV higher than the initial state and 0.035 eV higher than the physisorption state. The small energy 
difference of 0.035 eV also indicates that very large vacuum is needed if one wants to obtain the 
fully desorbed state in the vacuum. This calculated potential surface along the desorption MEP is 
very close to Hornekaer's theoretical result, 3 where the smoothly varying potential surface for the 
MEP from the physisorption state to the fully desorbed state also has been shown. In addition, if 
the hydrogen monomer climbs over the desorption barrier, it will escape from the graphene layer 
with high velocity after climbing down the potential barrier. In the diffusion process shown in Fig. 
Qlb), the hydrogen atom moves along a C-C bond at an average height of about 1.40 A. In the 
diffusion saddle-point state, the hydrogen atom is at the bridge-site with a height of 1 .20A. The 
diffusion potential barrier is 1.035 eV, which is somewhat smaller than the desorption potential 
barrier. 



The kinetic movement of the hydrogen monomer can result in the kinetic change of the mag- 
netic property of the hydrogenated graphene. Fig. [2] shows the electronic DOS of the initial, 
desorption saddle-point and diffusion saddle-point states, as well as the respective partial DOS of 
the hydrogen monomer. The initial state and the desorption saddle-point state both have a total 
magnetic moment of 1.0 \i B , while the diffusion saddle-point state is non-magnetic. The DOS of 
the initial state is consistent with the theoretical results by others. 20,45 The adsorption of hydrogen 
atom breaks an aromatic n bond in graphene, which results in one unsaturated C(p) orbital. This 
unsaturated C(p) orbital is responsible for the two narrow peaks (bands) placed around the Fermi 
level in the total DOS. The energy gap is 0.245 eV. As seen from the partial DOS of the H(ls) 
orbital, the contribution from the H(ls) orbital to these two bands is small. The unpaired electron 
occupying the lower band is responsible for the magnetic moment of 1.0 \Ib- The energies of the 
other spin-up bands shift downwards a little with respect to their spin-down counterparts due to 
the exchange-correlation interaction between the electrons in the magnetic C(p) orbital and the 
other electrons. In the desorption saddle-point state, the partial DOS of the H(ls) orbital is much 
larger than that in the initial state. This is because the H(ls) orbital is much less hybridized with 
graphene in the desorption saddle-point state than that in the initial state, and less electrons trans- 
fer from H to graphene. Thus in the desorption process, the contribution of the H(l s) orbital to the 
total magnetic moment of 1.0 \Ib will increase up to 100% when the H monomer is fully desorbed. 
The spin-up bands do not significantly shift with respect to the spin-down ones, which indicates 
the interaction between the magnetic H atom with graphene is weak. And the linear dispersion of 
the total DOS near the Fermi level is the same as that of grpahene, indicating that the aromatic n 
bond in graphene which is broken in the chemisorption state has been restored. In the diffusion 
saddle-point state, the linear dispersion part of the total DOS is not influenced by the presence of 
the hydrogen monomer on the bridge site and the Fermi level shifts upwards into the conduction 
bands. This reveals that the hydrogen monomer does not hybridize with graphene but just dopes 
graphene with itinerant electrons. This is the reason why the diffusion saddle-point state is non- 
magnetic. Thus in the diffusion process, the total magnetic moment decreases from 1.0 \Ib down 
to 0.0 \Xb when leaving from the chemisorption state to the diffusion saddle-point state. 

The displacement of the hydrogen monomer in each eigenvector of the vibrational dynamic 
matrix was represented by e(a» r ) (i is the index of the vibrational mode); the calculated spectra of 
the |e M (o>,)| 2 (M=H, D and T) for the initial, desorption saddle-point and diffusion saddle-point 
states are shown in Fig. [3j The modes with large values of |e(w,)| 2 are localized stretching (S) 
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and bending (B) modes. The isotope effect in the phonon spectra is obvious in that the vibrational 
frequency decreases with increasing monomer mass. The stretching modes of H at 2664 cm" 1 and 
of D at 1951 cm" 1 in the initial state are very close to the experimental measurements of 2650- 
2680 cm" 1 and 1950-1955 cm -1 .— & In the desorption saddle-point state, the stretching mode of 
each isotope monomer becomes imaginary, while in the diffusion saddle-point state, one bending 
mode of each isotope monomer becomes imaginary. The real stretching and bending modes in 
the diffusion saddle-point state have larger vibrational frequencies than the real bending modes in 
the desorption saddle-point state; this will result in a difference between the exponential prefactor 
of desorption and that of diffusion, as shown below. For light-mass atoms such as hydrogen 
isotopes, the vibrational zero-point energy corrections (AF vib (0)) are very important to the precise 
estimation of activation energies. The calculated AF vib (0)s for H, D and T monomers are listed 
in Tab. H The values, which range from 0.072 eV to 0.205 eV, indeed cannot be neglected if the 
kinetic properties of hydrogen isotope monomers are to be understood. The isotope effect in the 
AF v tb(0)s, which exhibits as that AF vib (0) decreases with increasing monomer mass, is a result of 
the isotope effect in the phonon spectra. Both of these isotope effects will result in the presence 
of an isotope effect in the desorption and diffusion frequencies of hydrogen monomers, as shown 
in the following. In addition, the AF vib (0) correction for the desorption potential barrier is larger 
than that for the diffusion potential barrier due to the disappearance of the highest-frequency real 
stretching modes and the significant lowering of the bending modes in the desorption saddle-point 
state with respect to those modes in the initial state. However, in the diffusion saddle-point state, 
only one lower-frequency real bending mode disappears and the real stretching and bending modes 
are less affected. The resonant in-band modes (modes with small values of |e(o»,)| 2 ) in these three 
states are slightly shifted and are nearly independent of the isotopic character of the adatom. The 
mode with a frequency around 630 cm" 1 observed in experiments is an example of the resonant 
in-band modes .-^^ 

The exponential prefactors for H, D and T monomers calculated from the quantum- 
mechanically modified hTST (v* m ) and its classical limit (v* cl , the Vineyard form—) are shown in 
Fig. HI The v*,s are much larger than the corresponding v*,„s in the plotted temperature range. 
Because differences in hydrogen monomer prefactors are mainly determined by the localized 
stretching and bending modes, and the contributions from the high-frequency localized modes are 
fully considered in the classical-limit hTST (Equ. H]), while they are much less considered in the 
quantum-mechanically modified hTST due to the weight factors of the bosonic phonon occupation 
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numbers (Equ. [3]). The quantum-mechanically modified and the classical-limit prefactors should 
be equal at very high temperatures. In the classical-limit hTST, all vibrational modes are regarded 
as completely thermo-activated. In the v*,s, the isotope effect is obvious: the prefactor decreases 
with increasing monomer mass, which is due to the isotope effect in the phonon spectra in Fig. 
[H as discussed in the previous paragraph. However, among the quantum-mechanically modified 
desorption prefactors, the relationship v* qm (T) > v* m (D) > v* m (H) holds and the isotope effect is 
reversed with respect to the isotope effect in the v* ; s. The values of the quantum-mechanically 
modified diffusion prefactors for H, D and T monomers are nearly the same over a very wide tem- 
perature range (0-1000K). This can be understood from Equ. [3] and the phonon spectra in Fig. [3j 
For convenience, the phonon spectra are divided into lower and higher regions below and above 
500 cm -1 , respectively. In the initial state, all the localized modes are in the higher region. The real 
localized bending modes in the desorption saddle-point state are significantly lowered such that 
they now fall within the lower region, while the real localized modes in the diffusion saddle-point 
state remain in the higher region. In the plotted temperature range, the v* m s are mainly determined 
by the vibrational modes in the lower region, due to the large bosonic phonon occupation numbers 
of these modes. In the phonon spectra of the initial and diffusion saddle-point state, the lower 
region contains only the resonant in-band modes, which are nearly independent of the isotopic 
character of the monomer, therefore, the v* s for H, D and T monomers in diffusion are nearly the 
same. However, in the phonon spectra of the desorption saddle-point state, the low-lying localized 
bending modes are significantly occupied by phonons at temperatures above 100 K; in this state, 
the localized modes for the T monomer are the lowest and should be occupied by most phonons 
and most readily thermo-activated, followed by the modes for D and H monomers in that order. 
This order of thermo-activation leads to the relative magnitude order of the v* m s for H, D and T 
monomers at temperatures above 100 K; however, this magnitude order will be progressively al- 
tered as the temperature increased to high values (not shown), finally achieving the same order as 



V cfi- 



The desorption frequencies (vd es ) and diffusion frequencies (ydiff) for the hydrogen-isotope 
monomers are shown in Fig. [51 where the v* qm is used as the prefactor. Both Vd es and Vdiff decrease 
with increasing monomer mass. This isotope effect is due to the isotope effect in the vibrational 
zero-point energy correction, as shown in Tab. [0 The v is exponentially dependent on the activa- 
tion energy and linearly dependent on the v* ; thus, the small reversed isotope effect in the v* s 
for the desorption (Fig. ffla)) does not compete with the normal isotope effect in the zero-point 



energy corrections. For each isotope, the desorption frequencies are larger than the corresponding 
diffusion frequencies because the desorption activation energy is less than the diffusion activation 
energy from Equ. [21 Additionally, the desorption prefactor is also larger than the diffusion prefac- 
tor. At room temperature (298 K), the v<#//S for H, D and T monomers are 1.6 x 10" 3 , 4.8 x 10~ 4 
and 3.0 x 10" 4 s~ : , respectively; the v des s are 1.4 x 10" 2 , 1.4 x 10~ 3 and 5.5 x 10" 4 s -1 , respectively. 
After a number of hydrogen monomers are deposited onto a graphene or graphite surface and 
if the system is annealed at a constant temperature T for a time interval to, the monomers are 
desorbed from the surface or diffuse away from the initial site. The kinetic properties of hydrogen 
monomers under this annealing condition were simulated by Equ. [6] and [91 with t taken to be 
10 minutes (600 s) in accordance with Hornekaer's experiment. 3 Fig. (6fa) shows the variations of 
the relative residual monomer number (n(t )/n(0)) and the diffusion radius (jdiff) with respect to 
T . The curves of n(to)/n(0) and r di ff present an isotope effect that a hydrogen monomer of lower 
mass is desorbed from the surface and diffuses on the surface more easily. During a 10-minute 
annealing, significant desorption does not occur in our simulation until T > 250 K. However, 
the monomers are nearly completely desorbed at T > 294 K for H, T > 312 K for D and To > 
320 K for T. Although the rj///(600s) increases rapidly with To, the monomers (H, D and T) on 
the graphene should be diffusionally immobile at any T because the diffusion radius for each 
hydrogen-isotope monomer is very short at temperatures below the complete-desorption tempera- 
ture, above which there will be no residual monomer left for measurement. This type of diffusional 
immobility matches experimental measurements obtained on D monomers. 318 The annealing pro- 
cess of D monomers on graphene at 298 K is visualized in Fig. [6](b) and (c) in order to compare 
with Hornekaer's experiment in Ref. [3]. In accordance with this experiment, the initial coverage 
(monomer number) is set to be 0.03% (n(0) = 114 in a lOOnm x lOOmn surface area) (Fig. [6tb)). 
After 10 minutes, only 33% of the monomers remain on graphene, a value that is comparable to 
the 20% observed in Hornekaer's experiment. n(to)/n(0) exponentially decreases with temperature, 
and the value of 20% corresponds to the annealing temperature of 301 K in this simulation. The 
discrepancy between simulation and experiment may arise from the fact that a small fraction of 
the monomers will be desorbed during the heating and cooling steps in the annealing process; the 
effect of STM measurement will also result in additional desorption of hydrogen monomers, which 
is not considered in this simulation; the interlayer interaction makes the hydrogen monomers eas- 
ier to escape from graphite surface in experiments. As mentioned in the methodology part, the 
interlayer interaction in bilayer graphene will lower the potential barriers by about 8 meV. This 
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is because the bilayer graphene is less flexible than the mono-layer graphene, 32 which makes the 
mono-layer graphene more reactive to bond with the hydrogen monomer. If a correction of -8 
meV is considered in the desorption activation energy, n(t )/n(0) equals 20% at 299 K. Thus, the 
interlayer interaction just results in temperature difference of couples of Kelvins, which indicates 
its weak effect on the kinetic properties of the hydrogen monomer and supports the validity of the 
mono-layer structural model used in this report. 

In some experimental measurements of kinetic properties, the temperature is increased at a 
constant rate (a). The kinetic properties of hydrogen monomers under this constant-rate heating 
condition were simulated by Equ. [[OlandQj], with a taken to be 1 .0 K/s. From Equ. \TT\ the relative 
desorption rate (Rd es (T, t), Fig. Ha)) is determined both by the desorption frequency (v des {T)) and 
the residual monomer number (n(t)). On heating, Vtu s (T) increases with temperature or elapsed 
time (Fig. |5]) and n(t) decreases (Fig. [7£b)); thus, there will be a desorption peak for each hydrogen- 
isotope monomer in the desorption spectrum of R des . The desorption peak for the H monomer is 
at 316 K, for the D monomer at 336 K and for the T monomer at 344 K. In the desorption spectra, 
an isotope effect can be seen that a hydrogen monomer of lower mass is desorbed from the surface 
more easily, which is the same as the isotope effect in the annealing process. The same as discussed 
in the previous paragraph, the weak interlayer interaction will only make the desorption spectra 
move to low temperatures by about 2 K. 

In addition, when the coverage of hydrogen atoms on graphene increases, hydrogen dimers 
will form and dominate on the surface . 2 ' 16117 In the desorption spectra under constant-rate heating 
(1.0 K/s), there are two peaks at 445 (490) K and 560 (580) K for H (D) dimers, which is totally 
different from the desorption spectra of H (D) monomer in Fig. [7J The kinetic properties of 
hydrogen dimers on graphene and the effect of dimer-dimer interaction have been investigated by 
us using nearly the same method, which will appear in a forthcoming paper. Generally, the isotope 
effects, localization nature of the vibrations of hydrogen dimers and the dependence of the kinetic 
properties on the localized modes are the same as those of the hydrogen monomer. This means the 
investigation on the kinetic properties of hydrogen monomer on graphene is instructive to other 
more complex situations. 
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IV. CONCLUSIONS 

The kinetic properties of hydrogen-isotope (H, D and T) monomers on graphene were simu- 
lated using a composite ab initio method consisting of density functional theory, density functional 
perturbation theory and harmonic transition state theory. The simulations were based on the po- 
tential barriers calculated from DFT and on the phonon spectra calculated from DFPT, which were 
used together to analyze calculated vibrational zero-point energy corrections, exponential prefac- 
tors (quantum-mechanically modified and classical-limit), jump frequencies and simulated kinetic 
properties. The chemisorption state and desorption saddle-point state both have a magnetic mo- 
ment of 1.0 |i B , while the diffusion saddle-point state is nonmagnetic. From the electronic DOS 
analysis, in the desorption process, the hybridization of the hydrogen monomer with graphene is 
reduced and the contribution of the H(ls) orbital to the total magnetic moment increases. And 
in the diffusion process, when leaving from the chemisorption state to the diffusion saddle-point 
state, this hybridization also decreases; the transferred electrons from H to graphene become more 
and more itinerant; the total magnetic moment decreases from 1.0 \i B down to 0.0 \i B - The spectra 
of the localized vibrational modes of hydrogen are hydrogen-mass-dependent, which results in the 
isotope effects in the kinetic properties (prefactor, activation energy, jump frequency, desorption 
rate and diffusion radius) of H, D and T monomers. The zero-point energy correction decreases 
with increasing monomer mass and is larger for desorption than for diffusion. This results in the 
isotope effect that a hydrogen monomer of lower mass is desorbed and diffuses more easily (with 
lower activation energies), and in that the desorption frequency of each hydrogen-isotope monomer 
is larger than its diffusion frequency. In the simulated 10-minute annealing of hydrogen-isotope 
monomers on graphene, the monomers are quite diffusionally immobile at temperatures lower than 
the complete-desorption temperature, and a large fraction (67%) of the D monomers on graphene 
will be desorbed at the annealing temperature of 298 K. These theoretical results are very close 
to those obtained by experimental observations of hydrogen adatoms on graphite surfaces. 318 The 
thermal desorption spectra (a = 1 .0 K/s) of the hydrogen-isotope monomers have also been simu- 
lated. The resulting desorption peak for H is at 316 K, for D is at 336 K and for T is at 344 K. The 
precise simulation of the kinetic properties of the hydrogen-isotope monomers on graphene de- 
scribed in this report can help to further precise studies of more complex hydrogenation processes. 
The precise theoretical predictions for radioactive tritium adatoms are especially useful because it 
is dangerous to conduct experiments on them in conventional laboratories. 
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TABLE I: Vibrational zero-point energy corrections (AF v u,(0)) in the activation energies (in eV) for the 
desorption (des.) and diffusion (dif.) of hydrogen-isotope (H, D and T) monomers. 



A^(0) 



A^(0) 



AF T , (0) 

via v ' 



des. 
dif. 



0.205 
0.114 



0.140 
0.084 



0.113 

0.072 



(a) 



saddle-point 

1.98 A 





FIG. 1 : (Color online) MEPs for the (a) desorption and (b) diffusion of a hydrogen monomer on a graphene 
layer. In the paths, the initial, saddle-point and final states are labeled with stars, and their structures are 
shown at the side of the figure. The final state in the desorption MEP is fixed to be the physisorption state. 
The yellow spheres represent carbon atoms and the smaller blue spheres represent hydrogen atoms. 
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FIG. 2: (Color online) The total electronic DOS and the partial DOS of the hydrogen atom of the (a) initial 
state, (b) desorption saddle-point state and (c) the diffusion saddle-point state. 
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FIG. 3: (Color online) The calculated spectra of the |e(a>,)| (i is the index of the vibrational mode) for 
the hydrogen-isotope monomers in the (a) initial, (b) desorption saddle-point and (c) diffusion saddle-point 
states. The imaginary vibrational modes in the saddle-point states are designated with negative frequencies. 
"S" and "B" represent the stretching and bending modes, respectively. The division of the lower region and 
higher region (see text) at 500 cm _1 is labeled by a gray dotted vertical line. 
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FIG. 4: (Color online) Variation of the exponential prefactors with respect to temperature for the (a) des- 
orption and (b) diffusion processes of the hydrogen-isotope monomers on graphene. The values from the 
quantum-mechanically modified hTST (v* ) are compared with their classical-limit values (v* d ). 
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FIG. 5: (Color online) Variation of the (a) desorption frequencies (vdes) an d (b) diffusion frequencies (v^//) 
of hydrogen-isotope monomers from the quantum-mechanically modified hTST with respect to the inverse 
of temperature. 
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FIG. 6: (Color online) (a) Variation of the diffusion radius and the relative residual monomer number with 
respect to the annealing temperature (?o = 600s). The three gray dashed vertical lines at the complete- 
desorption temperatures of H, D and T monomers (294, 312 and 320 K, respectively) are intended as 
a visualization tool, (b), (c) The visualization of the annealing process in (a) for D monomers at room 
temperature (298 K) for 10 minutes. The area of the graphene surface is lOOnm x lOOnm, and the initial 
coverage is set to be 0.03% («(0) = 1 14), according to the experiment in Ref. yfl. 
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FIG. 7: (Color online) (a) Simulated thermal desorption spectra for hydrogen-isotope monomers with the 
heating rate a = 1.0 K/s. (b) The residual monomer number during the thermal desorption. 
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